The Statistical Multifragmentation Model with Skyrme Effective Interactions 
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The Statistical Multifragmentation Model is modified to incorporate the Helmholtz free energies 
calculated in the finite temperature Thomas-Fermi approximation using Skyrme effective interac- 
tions. In this formulation, the density of the fragments at the freeze-out configuration corresponds 
to the equilibrium value obtained in the Thomas-Fermi approximation at the given temperature. 
The behavior of the nuclear caloric curve at constant volume is investigated in the micro-canonical 
ensemble and a plateau is observed for excitation energies between 8 and 10 MeV per nucleon. A 
kink in the caloric curve is found at the onset of this gas transition, indicating the existence of 
a small excitation energy region with negative heat capacity. In contrast to previous statistical 
calculations, this situation takes place even in this case in which the system is constrained to fixed 
volume. The observed phase transition takes place at approximately constant entropy. The charge 
distribution and other observables also turn out to be sensitive to the treatment employed in the 
calculation of the free energies and the fragments' volumes at finite temperature, specially at high 
excitation energies. The isotopic distribution is also affected by this treatment, which suggests that 
this prescription may help to obtain information on the nuclear equation of state. 

PACS numbers: 25.70.Pq, 24.60.-k, 31.15.bt 



I. INTRODUCTION 

Understanding the behavior of nuclear matter far from 
equilibrium, besides its intrinsic relevance to theoretical 
nuclear physics, is a subject of great interest to nuclear 
astrophysics, where the fate of supernovae or the proper- 
ties of neutron stars are appreciably influenced by the nu- 
clear equation of state (EOS) [J [US]. Thus, this area has 
been intensively investigated in different contexts during 
the last decades @, |, 0, 1 i, 0, H S El EH- Nuclear 
collisions, at energies starting at a few tens of MeV per 
nucleon, provide a suitable means to study hot and com- 
pressed nuclear matter [|, El El El El 51 EE El E3 • 
The determination of the nuclear caloric curve is of par- 
ticular interest as it allows one to infer on the existence 
of a liquid-gas phase transition in nuclear matter. Never- 
theless, owing to experimental difficulties, conflicting ob- 
servations have been made in different experimental anal- 
yses [H El El M M HI M M HE Ha HE SHI 11. 

Although there have been attempts to reconcile these re- 



sults [31(, this issue has not been settled. 



The properties of the disassembling system in central 
collisions, as well as the outcome of the reactions, have 
been found to be fairly sensitive to the EOS employed 
in the many theoretical studies using dynamical models 
that have been performed d El EH, Hill- However, m 
spite of their success in describing many features of the 



nuclear Multifragmentation process [32|, |33j, l34( , there has 
not been much effort to incorporate information based on 
the EOS in the main ingredients of statistical multifrag- 
mentation models. Yet, these models have recently been 
applied to investigate, for instance, the Isospin depen- 
dence of the nuclear energy at densities below the sat- 
uration value [35|, HI [37} ■ These calculations have sug- 
gested an appreciable reduction of the symmetry energy 
coefficient at low densities but other statistical calcula- 
tions [H, HI H(J indicate that surface corrections to the 
symmetry energy may also explain the behavior observed 
in those studies. Therefore, statistical treatments, which 
consistently include density effects, are most advisable 
for these studies. 

In this work, we modify the Statistical Multifragmen- 
tation Model (SMM) HI H3, HU and calculate some of 
its key ingredients from the finite temperature Thomas- 
Fermi approximation [3, HI 51 H3] using Skyrme effec- 
tive interactions. This version of the model is henceforth 
labeled SMM-TF. The internal Helmholtz free energies 
of the fragments arc calculated in a mean field approxi- 
mation, which is fairly sensitive to the Skyrme force used 
. This makes possible to investigate whether such sta- 
tistical treatments may provide information on the EOS. 
Furthermore, this approach allows to consistently take 
into account contributions to the free energy due to ex- 
citations in the continuum, in contrast to the traditional 
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SMM [49(. For consistency with the mean field treat- 
ment, the equilibrium density of the fragments at the 
freeze-out stage is also provided by the Thomas-Fermi 
calculations. Thus, in contrast with former SMM calcu- 
lations, fragments are allowed to be formed at densities 
below their saturation value. For a fixed freeze-out vol- 
ume, this leads to a systematic reduction of the free vol- 
ume, which directly affects the entropy of the fragment- 
ing system, the fragment's kinetic energies, and, also, 
the system's pressure. As a consequence, other proper- 
ties, such as the caloric curve and the multiplicities of the 
different fragment species produced, are also affected. 

We have organized the remainder of this work as fol- 
lows. In Sect.HIlwe discuss the modifications to the SMM 
and present the results obtained with this modified treat- 
ment in Sect. Hill Concluding remarks are drawn in Sect. 
IIV1 In Appendix [X] we provide a brief description of the 
Thomas-Fermi calculations employed in this work. 



charge, and energy conservation and thus the following 
constraints are imposed for each partition: 
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II. THEORETICAL FRAMEWORK 

In the SMM S El , it is assumed that a source, 
made up of Zq protons and Aq — Zq neutrons, is formed 
at the late stages of the reaction, with total excitation 
energy E* . This excited source then undergoes a si- 
multaneous statistical breakup. As the system expands, 
there is a fast exchange of particles among the different 
fragments until a freeze-out configuration is reached, at 
which time particle exchange ceases and the composition 
of the primary fragments is well defined. One then as- 
sumes that thermal equilibrium has been reached and 
calculates the properties of the possible fragmentation 
modes through the laws of equilibrium statistical me- 
chanics. A possible scenario consists in conjecturing that 
the freeze-out configuration is always attained when the 
system reaches a fixed pressure, i.e. the nuclear multi- 
fragmentation is an isobaric process. In this case, differ- 
ent statistical calculations predict a plateau in the caloric 
curve [12, HH, [H, [H, 0, HH ■ The situation is quali- 
tatively different if one assumes that, for a given source, 
the freeze-out configuration is reached at a fixed breakup 
volume V x . As found in many different calculations, a 
monotonic increase of the temperature with excitation 
energy takes place in this case |55l . [56l H3] • In what fol- 
lows we demonstrate that this is a consequence of the 
properties assumed for the fragments formed, and not of 
the fixed volume assumption. 

In this work we keep the breakup volume fixed for all 
fragmentation modes, and parametrize it through the ex- 
pression: 



In the above equations, -E| ource is the ground state energy 
of the source, Na,z denotes the multiplicity of fragments 
with mass and atomic numbers A and Z, respectively, 
Ba,z corresponds to the binding energy of the fragment, 
and e* A z (T) represents its excitation energy at temper- 
ature T. The Coulomb repulsion among the fragments 
is taken into account by the last two terms of the above 
equation which, together with the self ener gy c ontribu- 
tion included in Ba,z, give the Wigner-Seitz [58( approx- 
imation discussed in Ref. [4l|. The coefficient Ccoul is 
given in Ref. As discussed in Ref. (HI, the frag- 

ment's bindin g en ergy Ba,z is either taken from experi- 
mental values [60l | or is obtained from a careful extrapo- 
lation if empirical information is not available. The spin 
degeneracy factors, which enter in the calculation of the 
translational energy -Etransi are also taken from experi- 
mental data for A < 4. In the case of heavier fragments, 
this factor is neglected, i.e., it is set to unity for all nuclei. 

One should notice that the freeze-out temperature 
varies from one fragmentation mode / = {Na,z} to an- 
other, since it is determined from the energy conservation 
constraint of Eq. (|4|). Therefore, the average tempera- 
ture is calculated, as any other observable O, through 
the usual statistical averages: 
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where Sf denotes the entropy associated with the mode 
/. This entropy is calculated through the standard ther- 
modynamical relation 



V x = (l + x )V 



(1) 



where Vq denotes the volume of the system at normal 
density and x > is an input parameter. 

In the micro-canonical version of SMM, the sam- 
pled fragmentation modes [43| are consistent with mass, 
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where 



F = E-TS 
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is the Hclmholtz free energy. In the following, we write 
this quantity as 

F = Y1 N A , Z [-B AtZ + r A ,z(T) + flT(T)] + F Coul 

A,Z 

(8) 

where the contributions from the fragment's internal ex- 
citation (f A z ) and translational motion (f A a z s ) are ex- 
plicitly separated. The latter reads: 



In the above expression, At = y ^at * s ^ ne thermal 
wavelength, m n is the nucleon mass, gA,z is the spin 
degeneracy factor, and Vf denotes the free volume, i.e., 
it is the difference between V x and the volume occupied 
by all the fragments at freeze-out. The quantity -Fboui 
corresponds to the last two terms in Eq. 

Before we present the changes in the model associated 
with the Thomas-Fermi calculations, we briefly recall be- 
low the calculation of Hclmholtz free energy F in the 
SMM. 



A. The standard SMM 

In its original formulation [4l[ , the SMM assumes that 
the diluted nuclear system undergoes a prompt breakup 
and that the resulting pieces of matter collapse to nor- 
mal nuclear density, although being at temperature T. 
Therefore, the volume occupied by the fragments corre- 
sponds to Vq, so that 

Vf = X Vo ■ (10) 

The energy and entropy associated with the translational 
motion of the fragment are respectively given by: 

^tians ^trans , th trans / ~i ~i \ 

e A,Z — JA,Z + 1s A,Z — 2 1 ^ ii > 

and 



The values of the parameters in the above expression are 
en = 16.0 MeV, T c = 18.0 MeV and f3 = 18.0 MeV 
[49(. This expression is used for all nuclei with A > 
5. Lighter fragments are assumed to behave as point 
particles, except for the alpha particle, for which one 
retains the bulk contribution to the free energy in order 
to take its excited states into account. 

In Ref. [4!| , the calculation of f A z has been modified 
to include empirical information on the excited states of 
light nuclei. We label this version of the model as ISMM 
and it is used throughout this work. 

B. The SMM-TF 

The Thomas-Fermi approximation, briefly outlined in 
Appendix EJ allows one to calculate the internal free 
energy of the fragments f A z from Skyrme effective in- 
teractions. Equations (|A19IA26|) clearly show that f A z 
contains, besides those from the nuclear interaction tra- 
ditionally used in SMM, contributions associated with 
the Coulomb energy in addition to the ones appearing in 
[Eq. (| 13[) ] . The additional Coulomb contribution arises, 
in the present case, because the equilibrium density of 
the nucleus at temperature T does not correspond, in 
general, to its ground state value. This is illustrated in 
Fig. [TJ which shows the ratio between the average density 
(p) at a temperature T, and the corresponding ground 
state value (po), for several selected light and intermedi- 
ate mass nuclei. We define the sharp cutoff density (p) as 
that which gives the same root mean square radius as the 
actual nuclear density obtained in the Thomas-Fermi cal- 
culation. One observes that (p) decreases as one increases 
the temperature of the nucleus and that it quickly goes to 
zero as T approaches its limiting temperature, since the 
nuclear matter tends to move to the external border of 
the box due to Coulomb instabilities [I?], HH, IMI ■ In our 
SMM-TF calculations presented below, we only accept 
a fragmentation mode at temperature T if it is smaller 
than the limiting temperature of all the fragments of the 
partition. If this is not the case, the entire partition is 
discarded as not being physically possible and we sample 
another one. 

Thus, a fragment's volume at temperature T is defined 
as: 
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The internal free energy f A z has contributions from 
bulk and surface terms: 
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where V A z represents the ground state value. Since it 
is useful to have analytical formulae to use in practical 
SMM calculations, we performed a fit of (p A ' Z ) using the 
following expression 
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FIG. 1: (Color online) Ratio between the average equilibrium 
density of the nucleus at temperature T and the ground state 
value as a function of the temperature. For details, see the 
text. 



where {af' Z } are the fit parameters. This expression has 
proven to be accurate enough for numerical applications, 
as is illustrated in Fig. [TJ which shows a comparison be- 
tween Eq. (jT5j) (full lines) and the results obtained with 
the Thomas-Fermi calculation (circles). The fit was car- 
ried out using n — 6. The dashed lines emphasize the 
fact that {p A ' Z ) = (Pq' Z ) in the standard SMM. 

Instead of being given by Eq. (fTU)) , the free volume of 
a fragmentation mode now reads 



V f = (l+x)V - 
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For the values of x usually adopted in statistical calcu- 
lations (0 < x < 5), this expression shows that, for some 
partitions, there may be a temperature Ty for which 
V f < 0. Therefore, if Eq. g]) leads to T > T v , the 
partition is discarded as it is not a physically acceptable 
solution. 

From Eq. the entropy associated with the kinetic 
motion of the fragment (A, Z) becomes 
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One should notice that, besides the smaller free volume, 
the last term in the above expression does not appear 
in the earlier version of the SMM. Since dVf/dT < 0, 
the expression above gives a smaller contribution to the 
total entropy than Eq. fT2")) . Owing to this change in 
the entropy, the average kinetic energy of the fragment 
becomes: 
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FIG. 2: (Color online) Kinetic energy (full line) of a partic- 
ular partition of the 150 Nd nucleus into M = 29 fragments, 
for E* I A — 8 MeV and V x /Vo — 3, as a function of the tem- 
perature. For comparison, the standard average translational 
energy, (M — lOfT, is also displayed (dashed line). The dif- 
ference between the left and the right hand sides of Eq. Q, 
AE, is also shown (dashed-dotted line). 



which, for a given temperature, is also lower than the 
corresponding SMM value. As a matter of fact, if the 
second factor dominates the first one, for T > Tr-, where 
e^ a f s (T/f) = 0, it can even become negative. We also 
discard all partitions for which there is no solution of Eq. 
flU satisfying T <T K . This aspect is illustrated in Fig. [11 
which shows the total kinetic energy -Etrans as a function 
of T for the 150 Nd nucleus, with E*/A = 8 MeV, for 
a partition containing M = 29 fragments. The full line 
represents Straus, whereas the dashed line corresponds to 
the standard SMM formula. The factor, M — 1 is due to 
the fact that the center of mass motion is consistently 
removed in all the kinetic formulae, although it is not 
explicitly stated above. 

The observed drop of the kinetic energy may lead to 
nontrivial consequences. In the case of the 150 Nd nucleus 
and for E* /A < 7.0 MeV, the fragment multiplicity is 
relatively low. Therefore, in this lower excitation energy 
range the behavior of the kinetic energy does not lead to 
any qualitative changes arising from the energy conser- 
vation constraint. However, for higher excitation ener- 
gies, and consequently larger fragment multiplicities, the 
kinetic energy is comparable to the total energy of the 
system -Etotai- In this case, for a given value of -Etotai, 
there may be two values of T which are acceptable solu- 
tions to Eq. ((3]). This is also illustrated in Fig. [H which 
shows the difference AE between the left and the right 
hand sides of this expression. Since all the micro-states 
corresponding to the same total energy Etotai should be 
included, both solutions, in this case associated with tem- 
peratures T m 5.3 MeV and T « 6.2 MeV must be con- 
sidered. They contribute, however, with different statis- 
tical weights, due to the different number of states asso- 
ciated with each of these two solutions. 
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FIG. 3: (Color online) Temperature distribution for the 
breakup of the the 150 Nd nucleus at E*/A = 8 MeV and 
for V x /Vo = 3. The full line corresponds to all events while 
those in which there are two temperatures associated with 
Stotai are depicted by the dashed line. 

Based on this scenario, the determination of the freeze- 
out temperature from isotopic ratios [62l ] , where one tac- 
itly assumes that T is univocally determined from E*, 
should be carefully reexamined. To give a quantitative 
estimate of these effects, we show, in Fig. 02 the tem- 
perature distribution for the fragmentation of the 150 Nd 
nucleus, at E*/A = 8 MeV and V x = 3V . The full line 
in this picture shows the results when all the partitions 
are considered whereas the dashed line represents only 
those which lead to two different temperatures. The cases 
where there are two temperature solutions correspond to 
43% of the events and account for 76% of the total statis- 
tical weight. These numbers are drastically changed at 
lower excitation energies where, for instance, one finds, 
at E*/A = 6 MeV, 0.03% and 0.09%, respectively. In 
spite of the great importance of these solutions at high 
excitation energy, the temperature distribution does not 
exhibit two clear dominant peaks, separated by a gap, as 
it could be expected from Fig. [2j This is because the nu- 
merical value of the two solutions vary from one partition 
to the other and the expected signature is thus blurred. 

We have also fitted the internal free energies of the 
nuclei through a simple analytical formula: 

TCI 

flz = -T 2 Y, b t' ZT ^ ( 19 ) 

i=0 

where {bf' Z } are the fit coefficients. The results are de- 
picted in Fig. 0] by the full lines, whereas the Thomas- 
Fermi calculations are represented by the circles. As 
in the previous case, an excellent agreement is obtained 
with a small number of parameters (to = 5). The free 
energies used in the ISMM are also shown in this pic- 
ture and are represented by the dashed lines. One sees 
that there are noticeable differences at low temperatures, 
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FIG. 4: (Color online) Internal free energy of selected nuclei 
as a function of the temperature. For details, see the text. 

in the case of the lighter nuclei. Particularly, many more 
states are suppressed in the ISMM than in the SMM-TF, 
which suggests that the latter should predict larger frag- 
ment multiplicities than the former. This is due to the 
empirical information on excited states which are taken 
into account in the ISMM [49[. In the case of heavier 
nuclei, the differences are more important at higher tem- 
peratures where the ISMM has more contributions from 
states in the continuum than the SMM-TF. However, the 
determination of the free energy at high temperatures in 
the ISMM is not as reliable as in the Thomas-Fermi ap- 
proximation in the sense that the numerical values of the 
parameters eo, T c , and /?o, used in actual calculations, 
are not obtained from a fundamental theory. They cor- 
respond to average values [4l|, [52] which, sometimes, are 
slightly changed by different authors (42|, 5^, [63| . 

From the above parametrization to f* A z , the en- 
tropy and excitation energy associated with the fragment 
(A, Z) read: 



s*a,z = 2T^bf' Z T l +T 2 ^ibf' Z T 1 - 1 (20) 

4=0 i=X 

and 

m m 

<A.z = t2 J2 h i' ZTi + TZ H 'V''/' : • (21) 

4=0 4=1 

The free energies and equilibrium volumes are calcu- 
lated using the above expressions for alpha particles and 
all nuclei with A > 5. 



III. RESULTS AND DISCUSSION 

The SMM-TF model described in the previous section 
is now applied to study the breakup of the 150 Nd nu- 
cleus at fixed freeze-out density. We use V x /Vq = 3 
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FIG. 5: (Color online) Caloric curve associated with the 
breakup of the 150 Nd nucleus. The ISMM calculation of Ref. 
[4gj ] (triangles) and the SMM-TF calculation presented in this 
work (circles) are seen to differ for E* /A > 8 MeV. For refer- 
ence the excitation energy of the compound nucleus calculated 
within the Fermi gas model (full line) , the classical gas model 
(dashed line) as well as the Thomas-Fermi approach (dotted 
line) have also been presented. For further details, see the 
text. 



in all calculations below. The caloric curve of the sys- 
tem is displayed in Fig. [5] Besides the SMM-TF (cir- 
cles) and the ISMM (triangles) results, the Thomas- 
Fermi calculations for the 150 Nd nucleus are also shown 
(dotted line), as well as the Fermi gas (full line) and 
the Maxwcll-Boltzmann (dashed line) expressions. For 
E*/A < 8.0 MeV, both SMM calculations agree fairly 
well on the prediction of the breakup temperatures. How- 
ever, a kink in the caloric curve is observed at this 
point, in the case of the SMM-TF, indicating that the 
heat capacity of the system is negative within a small 
excitation energy range around this value. Negative 
heat capacities have been predicted by many calculations 
and have been stron g l y d ebated in the recent literature 
[H [H [EH, [H, HI HI M, HE HI M, HI . However, this 
feature is normally observed at the onset of the multi- 
fragment emission, i.e. at the beginning of the liquid-gas 
phase transition [42l . |53| , whereas it appears much later 
in the present calculation. 

In order to understand the qualitative differences be- 
tween the two SMM approaches, we show, in Fig. [51 
the multiplicity of light particles Ni p (all particles with 
A < 4, except for alpha particles), the alpha particle N a 
and the Intermediate Mass Fragment (IMF, 3 < Z < 15) 
Njmf multiplicities, as well as the total fragment mul- 
tiplicity -ZVtotai as a function of the excitation energy. 
It is important to notice that neutrons are included in 
Nip and iVtotai- One sees that there is a clear disagree- 
ment between the two SMM calculations in the predic- 
tion of the alpha particle multiplicity. This is due to the 
construction of the internal free energies in the ISMM 
[4t|, which considers empirical low energy discrete states. 
Since the first excited state of the alpha particle is around 



FIG. 6: (Color online) Average multiplicity of light particles, 
alphas, IMF's and the total fragment multiplicity, as a func- 
tion of the excitation energy. For details, see the text. 



20 MeV, this strongly increases the free energy at low 
temperatures within the ISMM calculation, in contrast 
to the Thomas-Fermi model calculations. In the case of 
the other multiplicities, the agreement between the two 
model calculations is fairly good for excitation energies 
up to E* I A w 8 MeV. The small discrepancy between 
ZVtotai in the two calculations can be attributed to the 
differences in the alpha multiplicities. All multiplicities 
rise smoothly up to approximately this excitation energy. 
Then, at E*/A w 8 MeV, in the SMM-TF calculations, 
N a and Nimf reach a maximum and decrease from there 
on. This behavior is not observed in the case of the ISMM 
because it takes place beyond the energy range consid- 
ered in the figure. Another feature also observed in this 
picture is the sudden change in the slope of the iV tota i and 
Nip multiplicities calculated using the SMM-TF model, 
which also takes place at the excitation energy mentioned 
above, and which is not seen in the ISMM results. 

Although the Hclmholtz free energies of the fragments 
are somewhat different in both calculations, the differ- 
ences are not large enough to quantitatively explain this 
peculiar behavior, as illustrated in Fig. 2J Therefore, 
the differences in the multiplicities calculated within the 
ISMM and SMM-TF models must be associated with the 
behavior of the kinetic terms, due to changes in the free 
volume in the SMM-TF calculations. 

To examine this aspect more closely, we show, in Fig. 
O the energy dependence of (Vf). It confirms the expec- 
tation that (Vf) should decrease as E* increases, ow- 
ing to the expansion of the fragments' volumes at fi- 
nite temperature. However, it reaches a minimum at 
E*/A w 8.0 MeV and rises from this point on. The 
logarithmic volume term of the entropy [Eq. (TIT)) ] dis- 
favors partitions with small free volumes. Furthermore, 
the last term in Eq. (TIT)) also gets larger as T increases 
since, besides being explicitly proportional to T, the fac- 
tor | -4^ | grows faster at high temperatures, as it can 
be inferred from the behavior of the densities shown in 
Fig. [T] Therefore, the system favors the emission of very 



light particles, Ni p , which cannot become excited in our 
treatment, in order to minimize the reduction of Vf. Nev- 
ertheless, this preference is closely related to the energy 
conservation constraint given by Eq. It is only when 
the excitation energy becomes sufficiently high that there 
is enough energy for the system to produce a significant 
number of very light particles. The inset in Fig. [7] shows 
the entropy per nucleon predicted by the two SMM treat- 
ments. It reveals that, while in the ISMM case it rises 
steadily, the entropy saturates, and even decreases in the 
SMM-TF model for 8.0 < E*/A < 11.0 MeV. The large 
emission of particles which have no internal degrees of 
freedom prevents the entropy from falling off from this 
point on, since they lead to larger (smaller absolute 
values) by increasing Vf, as they do not expand. One 
should notice that the reduction of the complex fragment 
multiplicities does not mean that the limiting tempera- 
ture of the fragments in the different partitions has been 
reached. In fact, the breakup temperatures obtained in 
the present calculations are much lower than the limiting 
temperatures of most nuclei, except for the very asym- 
metric ones, as may be seen in the examples given in 
Fig. [T]and in Refs. [4]], El]. This effect on the fragments 
produced should appear at much higher excitation ener- 
gies, as those fragments have excitation energies much 
smaller than the original nucleus, since an appreciable 
amount of energy is used in the breakup of the system. 
Therefore, the back bending of the caloric curve and the 
small plateau observed in Fig. [5] arc strongly ruled by 
the changes in the free volume. As a consequence of this 
fact, the phase transition at high excitation energy takes 
place at approximately constant entropy. 

This observation is also corroborated by the charge 
distributions shown in Fig. [5] for four different excitation 
energies: E*/A = 5, 6, 7 and 8 MeV. It shows that the 
multiplicity of heavy fragments is strongly reduced in the 
SMM-TF calculations as the excitation energy increases, 
although they are not completely ruled out of the pos- 
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sible fragmentation modes. In particular, the SMM-TF 
model systematically gives much lighter fragments than 
the ISMM, for the reasons just discussed. 

Even though the fragments are not directly affected 
by their limiting temperatures at the excitation energies 
we consider, the reduction of the entropy associated with 
the volume affects the fragment species in different ways. 
Indeed, since the proton rich nuclei tend to be more un- 
stable, they arc hindered due to these dilatation effects 
more strongly than the other isotopes. Owing to their 
larger volumes at a given temperature T, partitions con- 
taining proton rich fragments have smaller entropies than 
the others. Therefore, one should expect to observe a re- 
duction in the yields of these fragments. This qualitative 
reasoning is confirmed by the results presented in Fig. [51 
which displays the isotopic distribution of some selected 
light fragments, produced at E* /A = 6.0 MeV. One sees 
that, even though both SMM models make similar pre- 
dictions for many observables at this excitation energy, 
the role played by the free volume effects just discussed is 
non-negligible. Since the limiting temperatures, as well 
as the equilibrium density at temperature T, is sensitive 
to the effective interaction used [47], El] , these findings 
suggest that careful comparisons with experimental data 
may provide valuable information on the EOS. 



IV. CONCLUDING REMARKS 

We have modified the SMM to incorporate the 
Helmholtz free energies and equilibrium densities of nu- 
clei at finite temperature from the results obtained with 
the Thomas-Fermi approximation using Skyrme effective 
interactions. Owing to the reduction of the fragments' 
translational energy at finite temperature, the model pre- 
dicts the existence of two temperatures associated with 
the same total energy. This feature is directly associated 
with the reduction of the free volume due to the expan- 
sion of the fragments' volumes. If this statistical treat- 
ment proves to be more appropriate to describe the nu- 
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clear multifragmcntation process than its standard ver- 
sion, the determination of the isotopic temperatures, at 
high excitation energies, should be carefully reexamined, 
since one tacitly assumes a univocal relationship between 
the temperature and the excitation energy in the deriva- 
tion that leads to the corresponding formulae [62| . 

The thermal dilatation of the fragments' volumes also 
has important consequences on the fragmentation modes. 
For excitation energies larger than approximately 8 MeV 
per nucleon, it favors enhanced emission of particles 
which have no internal degrees of freedom (very light nu- 
clei, protons and neutrons), leading to the onset of a gas 
transition at excitation energies around this value. The 
existence of a small kink in the caloric curve, as well as a 
plateau, for a system at constant volume is qualitatively 
different from the results obtained in previous SMM cal- 
culations where these features were observed only at (or 
at least at nearly) constant pressure [5{|. 

Since many-particle multiplicities, such as those asso- 
ciated with the IMF's and the light particles, are very 
different in both statistical treatments for excitation en- 
ergies larger than 8 MeV per nucleon, we believe that 
careful comparisons with experimental data may help to 
establish which treatment is better suited for describing 
the multifragmcnt emission. Furthermore, since the iso- 
topic distribution turns out to be sensitive to the treat- 
ment even at lower excitation energies, this suggests that 
one may obtain important information on the EOS by 
using different Skyrme effective interactions in the SMM- 
TF calculations. Particularly, this modified SMM model 
is appropriate to investigate the density dependence of 
the symmetry energy recently discussed [3^, l36l . [37], E3| ■ 
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APPENDIX A: THE FINITE TEMPERATURE 
THOMAS-FERMI APPROXIMATION 

The Thomas-Fermi approximation to nuclear systems 
is thoroughly discussed in Refs. 0, [H, [H, Thus, 
we review its essential features below in order to give a 
full account of all calculations presented in this work. 

The equilibrium configuration of a nucleus at temper- 
ature T is found by minimizing the thermodynamical po- 
tential f2 with respect to the number density p a (a = p,n 
for protons or neutrons): 

n = F[p]-Y^ I d 3 ffi aPa (Al) 

a J 

where the Hclmholtz free energy is given by 



Wnucl + ^Co 



d A r 



(A2) 



In the above expression, S a denotes the entropy density 
associated with the species a, fi a is the corresponding 
chemical potential, 7i n uci is the nuclear energy density of 
the system, and the Coulomb term reads: 



Wcoui = yP P (r) 



4 V 7T 



1/3 



Pt /3 (rl 



(A3) 



The second term above corresponds to an approximation 
to the exchange contribution to the Coulomb energy 

m - n 

The expression for 7i nuc i given in Ref. [44J may be 
rewritten as 



Wnucl — Hq + Ht + Wgrad + Hj 



where 



Ho — 



to 
2 

*3 

12' 



x 



:i + y)p 2 -(-3 + i)(p J 2 l + ^) 



7~t T = T r 



(A4) 



(A5) 



2m! 



2m! 



(A6) 
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where 



rad 



l 

64 
1 

64 



9h - 5t 2 (l + -x 2 ) 
5 



VP 



(A7) 



[3ti(l + 2xi) + * a (l + 2x 2 )] ( VPn - VP P 



Hj = -W 



J • VP + Jn- VPn + Jp ■ VPp 



(A8) 



the total density is denoted by p = p n + p p and J = 
Jp + J n is the spin-orbit density. The kinetic factor r a is 
given by 



1 

2^2 



2< 



5/2 



r 5/2 / 3/2 (y a ) 



(A9) 



where 



2m* 



_9_ 

<9r Q 

__. 

2m 
1 r 
4 



'^nucl 

i 



+ - - xi) + 3t 2 (l + X 2 )] pa 

p a > (A10) 

and p a i = p p (p n ) if ct = n (p). The Fermi-Dirac integral 



l + y)+t 2 (l + y) 



l n /2(y) 



dx 



r n/2 



1 + CXp(.T — J/) 



(All) 



is efficiently calculated using the formulae given in Ref. 
fn\ . where one also finds approximations to the inverse 
function y(I n / 2 ). The latter is determined from the num- 
ber density 



1 /2m* 



3/2 



T 3/2 hMy a ) 



2^ 2 V h 2 

The entropy density S a can then be easily calculated 



(A12) 



5 h 2 



3 2ml T 



PaVa 



(A13) 



The parameter set {xi,ti,<j,Wo}, i = 0,1,2,3, for the 
Skyrme interaction used in this work, SKM, is listed in 
Ref. [4~jj . Since we stay in the zero-th order approxima- 
tion in h, J a = and then Tij does not contribute to 

H nucl M- 

Following Suraud and Vautherin [46|, S3 , the equilib- 
rium configuration is found by iterating the densities at 
the k-th step according to 



,(k+x) 



Pi k) 



Pa 



(A14) 



,(*0 



N p = Z, N n = A 



dV£«(r>«(r), 



(A15) 



Z. and 



ST 

Spa 



(A16) 



The parameter A is chosen to be small enough in order 
to ensure that the first order approximation given by Eq. 
(IA14I) remains valid. 

In our numerical implementation, we have assumed 
spherical symmetry, and discretized the space using a 
mesh spacing AR = 0.1 fm, which suffices for our pur- 
poses. As suggested in Refs. 0, [47[, the second term of 
Eq. (|A7[) is neglected since it is small and may lead to nu- 
merical instabilities. Similarly to the treatment adopted 
in Ref. [46| , the gradient density terms are calculated at 
the mesh point r i+l / 2 = (i+ 1/2) AR, using [72| 



T^p{r i+1/2 ) 



p{rj) - p(rj_i) 
2(Ai?/2) 



0[{AR) 2 



(A17) 



which turned out to be numerically stable. 

Due to the important contributions associated with 
unbound states at high temperatures, the above treat- 
ment is not accurate for T > 4 MeV, as pointed out 
by Bonche, Levit, and Vautherin [6l[. Therefore, these 
authors have proposed a method to extend the Hartree- 
Fock calculations to higher temperatures. As they have 
noticed, there are two solutions of the Hartree-Fock equa- 
tions for a given chemical potential. One of them cor- 
responds to a nucleus in equilibrium with its evaporated 
particles whereas the other is associated with the nucleon 
gas. Thus, in their formalism, the properties of the hot 
nucleus is obtained by subtracting the thermodynamical 
potential associated with an introduced nucleon gas 
from that corresponding to the nucleus in equilibrium 
with its evaporated gas Ung- Except for the Coulomb 
energy, there is no interaction between the gas and the 
nucleus-gas system. 

This approach has been successfully applied by these 
authors jH, [6l| and has been adapted to the finite tem- 
perature Thomas-Fermi approximation by Suraud [47| . 
More precisely, the thermodynamical potential associ- 
ated with the nucleus is given by 



E, 



Coul 



(A18) 



One should notice that, by construction, Qng an d do 
not contain any Coulomb contribution. More specifically, 
one defines the subtracted free energy 
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since p NG and p G are constrained by 



■sub 



rl nucl — rl nucl — J V a a 



d 3 i 



K 



sub j3 



Coul 



<Tr 



(A19) 



<V't 



pt NG) (r) 



p^ G \r) 



(A24) 



where the subtracted Coulomb energy density, in the last 
term of this expression, reads 



ft 



Coul 



z j \ r — r 



4 V 7T 



1/3 



« G ) 4/3 - (P P G ) 4/3 



and the subtracted density p p : 



(A20) 



One then starts with a reasonable guess for p^ G and 
which can be a Woods-Saxon density for the for- 
mer and a small constant value for the latter (subject 
to the condition p a > 0), obeying the constraint given 
by the above expression, and apply the iteration scheme 
just described. Ideally, convergence is reached when 

Ba (0 — vanishes, so that pi fe ' 7 ^ becomes station- 
ary. In practice, one can monitor the quantity 



P G 



P P (rl = p^(f)-p G (f) 



(A21) 



is the quantity that enters in the direct part of the 
Coulomb energy. 

The iteration scheme given by Eq. (I A14|) remains un- 
changed if one rewrites B)£' as 



sub 



5p^ ] 



(A22) 



where the super- index (k, 7) denotes the quantity associ- 
ated with the gas (7 = G) or the nucleus-gas (7 = NG) 
at the &-th stage of the iteration. The positive sign is 
associated with the NG solution whereas the negative 
sign is used in the other case. The proton and neutron 
chemical potentials are given by an expression similar to 
Eq. (IAT51) 



- Bi k - G \f)pi k - G Hr)}, (A23) 



A£ = 



d 3 r 



+ (^ G )(r)- M W)%i fe ' G )}(A25) 



and stop the iteration when the established tolerance is 
reached. The Hclmholtz free energy of the nucleus can 
then be easily calculated through Eq. (|A19[) . so that the 
internal free energy of the nucleus is 



/* Z (T) = jP sub (T) - T snb {T = 0) 



(A26) 



We have used the approximation just described in this 
Appendix to calculate j* A z for all the fragments entering 
in the SMM, with A > 5 (and alpha particles) from T = 
MeV up to the limiting temperature [13, El, [HI in steps 
of 0.1 MeV. 
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